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O ■ ABSTRACT 

o 

fH I We analyze the ASCA spectrum of the Cen X-3 X-ray binary system in echpse using 

Qh' atomic models appropriate to recombination-dominated level population kinetics in an 

O ' overionized plasma. In order to estimate the wind characteristics, we first fit the eclipse 

-4—' ' spectrum to a single-zone photoionized plasma model. We then fit spectra from a range 

C^ ■ of orbital phases using global models of photoionized winds from the companion star and 

the accretion disk that account for the continuous distribution of density and ionization 
k> , state. We find that the spectrum can be reproduced by a density distribution of the 

H I form derived by Castor, Abbott, & Klein (1975) for radiation-driven winds with Af /t>oo 

consistent with values for isolated stars of the same stellar type. This is surprising 
because the neutron star is very luminous (-^ 10'^^ erg s^^) and the X-rays from the 
neutron star should ionize the wind and destroy the ions that provide the opacity for 
the radiation-driven wind. Using the same functional form for the density profile, we 
also fit the spectrum to a spherically symmetric wind centered on the neutron star, a 
configuration chosen to represent a disk wind. We argue that the relatively modest 
orbital variation of the discrete spectrum rules out a disk wind hypothesis. 



1. Introduction 

Since the early observations of the High Mass X-ray Binary (HMXB) Cen X-3, it has been 
known that the system exhibits a residual X-ray flux in eclipse, indicating that the X-rays are 



scattered or otherwise reprocessed, as might occur in a wind (Schreier et al. 1972). A similar 
phenomenon has been observed in other HMXBs — for example, Becker et al. (1978) observed a 
residual eclipse flux from Vela X-1. It has been realized since the far ultraviolet became accessible 
with rocket and satellite-borne instruments that isolated hot stars expel material at rates of order 
10~® Mq yr~^ (Morton 1967), and so the existence of winds in HMXBs might have been expected. It 
was shown that in isolated hot stars, strong winds are driven by transfer of the outward momentum 
of the ultraviolet stellar radiation to the matter through line absorption (Lucy & Solomon 1970; 
Castor et al. 1975). 

However, it was pointed out by Hatchett & McCray (1977) that in an X-ray binary. X-rays from 
the compact object would ionize a part of the wind around the X-ray source to such a high degree 
that the ions with line transitions which enable radiative driving would not be present. Presumably, 
if an X-ray source were luminous enough, it could completely shut off the radiation driven wind 
on the X-ray illuminated side of the companion star. Calculations by MacGregor & Vitello (1982) 
showed that radiative driving is disabled for X-ray luminosities larger than ~ 5 x lO^^ergs"^. 
As HMXBs typically have luminosities in the range 10^^-10^^ erg s~^, this would indicate that 
radiative driving could not function at all in X-ray binaries. However, this calculation assumed 
that the wind is optically thin to the ionizing X-rays. The effects of the wind's opacity to X- 
rays has been explored by Masai (1984) and by Stevens (1991). These calculations considered the 
formation of a He^"'"/He^ Stromgren surface in the wind. On the He"*" side of this boundary, all 
elements have lower ionization, and more of the ions necessary for radiative driving are present. 
Masai (1984) concluded that most wind- fed X-ray binaries must contain Stromgren boundaries and 
hence radiative winds. This conclusion is limited to wind- fed systems because of the coupling of 
the wind parameters to the accretion luminosity in these systems. An increase in the luminosity 
of the X-ray source tends to increase the volume of the He^"*" region, however, to increase the 
luminosity it is necessary to increase the mass loss rate of the companion star or decrease the 
wind velocity, both of which tend to decrease the volume of the He^"*" region. The investigation 
of Stevens (1991) was not limited to wind-fed systems. It showed that X-ray luminosities as large 
as ~ lO^^ergs"^ diminished the wind velocity and mass loss rate but did not completely shut 
off the wind. However, for luminosities larger than 4 x lO^^ergs"^ Stevens (1991) was "unable 
to to find dynamical solutions" for the wind. Ho &: Arons (1987) explored the behavior of high 
luminosity (exclusively) wind-fed X-ray binaries. However, they did not explicitly consider optical 
depth effects and their condition for X-ray ionization turning off the wind is questionable, even for 
the optically thin case (c.f. Stevens 1991). We regard the nature of winds in X-ray binaries outside 
of the X-ray shadow of the companion as unresolved, and especially problematic in the case of high 
luminosity (Lx ~ lO^^ergs"^) systems. 

A wind driving mechanism which does not depend on the presence of ions with UV resonance 
lines is thermal pressure due to X-ray heating of the exposed face of the companion star (also 
referred to as evaporative winds). This mechanism was invoked by Basko & Sunyaev (1973) and by 
Arons (1973) to explain the mass transfer in Her X-1, though Alme & Wilson (1974) and McCray 



& Hatchett (1975) found that a thermal wind alone could not power the X-ray source. Day &: 
Stevens (1993) showed that an X-ray excited wind could account for the mass transfer in Cen X-3, 
as well as explain the extended eclipse transitions seen in the source. At least one simulation of 
the disk-fed high-mass X-ray binary LMC X-4, in which the effects of X-ray heating were included, 
showed that a stronger wind was driven from the accretion disk, and that the structure of the wind 
was dominated by the disk wind (Owen & Blondin 1997). 

Before the launch of ASCA in 1993 the energy resolution of most cosmic X-ray detectors was 
rather poor (AE/E --^ 10-20%) and therefore, these instruments were unable to detect the X-ray 
emission lines which are a signature of optically thin, highly ionized gas. Nor is there any other 
waveband in which discrete emission from very highly ionized gas could be detected. Most obser- 
vational X-ray examinations of winds in HMXBs focused on absorption in low-ionization material, 
iron fluorescence, which is produced primarily in low-ionization material, and the Compton scat- 
tering continuum which is insensitive to the ionization state. (Kallman & White 1982; Nagase 
et al. 1986; Sato et al. 1986; Clark et al. 1988; Haberl et al. 1989; Woo et al. 1995). The lack of a 
regular low-ionization wind on the X-ray illuminated side of HMXBs, and therefore the presence of 
a high-ionization wind there, has been inferred from decreases in equivalent widths and velocities 
of the P Cygni profiles of ultraviolet resonance lines away from X-ray eclipse (Dupree et al. 1980; 
van der Klis et al. 1982; Hammer schlag-Hensberge et al. 1984; Vrtilek et al. 1997). However, direct 
study of the highly ionized material requires X-ray spectroscopy. The Solid-state Imaging Spec- 
trometers (SIS, Gendreau 1995) on board ASCA have energy resolution of a few per cent which 
allows identification and study of many previously undetectable X-ray spectral features. With the 
SIS detectors, recombination and fluorescence emission features were seen from several ions from 
the HMXBs Vela X-1 (Nagase et al. 1994), Cen X-3 (Ebisawa et al. 1996), and Cyg X-3 (Kitamoto 
et al. 1994; Liedahl & Paerels 1996; Kawashima & Kitamoto 1996). 

When a compact X-ray source is occulted by its companion, the emission spectrum from an 
extended wind can be studied without the confusion from the more intense, generally featureless, 
spectrum of X-rays from the neutron star. Sako et al. (1999) studied the eclipse spectrum of the low 
luminosity HMXB Vela X-1 obtained by ASCA and estimated the rate of mass loss. Though a highly 
ionized wind exists in Vela X-1, they found that that most of the mass is inside dense clumps, which 
are not highly ionized. This allows for the possibility that radiation imparts its outward momentum 
to the clumps which then drag the hot, diffuse wind outward. Presumably, these clumps could be 
destroyed (or inhibited from forming) by a more luminous X-ray source such as Cen X-3 or SMC X- 
1. In fact, Ebisawa et al. (1996) showed that the equivalent width of the 6.4 keV iron fluorescence 
line in Cen X-3 was nearly constant with orbital phase, indicating that most of the low-ionization 
material in the system is located near the neutron star, and that very little low-ionization material 
is found in the extended wind. This is also confirmed by pulsations in the 6.4 keV line (Day et al. 
1993; Audley 1998). Wojdowski, Clark, & Kallman (2000) used hydrodynamic simulations of the 
wind in the most luminous persistent HMXB SMC X-1 by Blondin & Woo (1995) to predict the 
X-ray eclipse spectrum of that system and compared it to a spectrum obtained with ASCA. In the 



simulation, a tenuous, very highly ionized wind formed on the X-ray illuminated side and a denser 
wind developed on the X-ray shadowed side. However, dense finger-like structures protruding from 
the shadowed side of the companion were swept into the X-ray illuminated region by the Coriolis 
force. The calculations showed that the reprocessed radiation from the tenuous gas was dominated 
by Compton scattering and the denser gas emitted copious amounts of recombination radiation. 
The observed spectrum was nearly featureless however, and Wojdowski et al. (2000) concluded that 
the dense fingers that appeared in the simulation could not be present in the wind of SMC X-1. 

Cen X-3 is one the most luminous, persistent known HMXBs in the Galaxy. It consists of a 
4.8 second pulsar in a 2.08 day eclipsing orbit (Schreier et al. 1972) with its 06-8 III type (Conti 
1978; Hutchings et al. 1979) companion V779 Cen. The high X-ray luminosity of Cen X-3 makes 
it an excellent candidate for the study of X-ray photoionized winds. It was observed by ASCA 
over approximately half an orbit, which included an eclipse. Ebisawa et al. (1996) found several 
emission lines in this data set, which were mostly from hydrogen-like ions. From the intensities of 
these lines, they made estimates of the scale and physical conditions of the wind. 

We re-analyze the data set of Ebisawa et al. (1996) using the observed spectra to test physically 
motivated wind models with the goal of providing constraints on the wind driving mechanism. In 
§ 2 we describe our reduction of the data. In § 3, we calculate emission spectra for photoionized 
plasmas using a list of ~3000 lines and emission features, fit the observed spectra using single zone 
emission spectra and then use the results to estimate wind parameters. In § 4 we calculate spectra 
using explicit parameterized models of the wind density distribution, and fit the observed spectra to 
determine wind parameters more accurately. We test two explicit global wind models: 1) a stellar 
wind from the companion with the velocity profile of a radiatively driven wind and 2) motivated by 
the accretion disk wind in the simulation of Owen & Blondin (1997), a wind with the same velocity 
profile but centered on the neutron star. In § 5 we justify an assumption that the wind is optically 
thin to X-rays which is used in previous sections. In § 6, we discuss the implications of our results. 



2. Data Reduction 

We obtained the screened REV2 ASCA event data from the 1993 June 24-25 observation of 
Cen X-3 across an eclipse from the HEASARC archive. All manipulation of the data was done 
with FTOOLS v4.2 (The Ftools Group 1998) and all of the programs mentioned in this section 
are from that package. We divided the data into the same four time segments as Ebisawa et al. 
(1996), corresponding to phase ranges -0.31 to -0.29, -0.23 to -0.08, -0.08 to 0.13, and 0.14 to 
0.20. The data was taken in a mix of FAINT and BRIGHT modes. For the data taken in FAINT 
mode, we used the files which had been converted to BRIGHT mode on the ground in order to 
have homogeneous data. For each of the time segments, we extracted all counts from inside a circle 
of radius 191". We extracted a background spectrum from an annulus of inner radius 191" and 
outer radius 382". During the observation, the center of the image was placed so that the source 
counts were distributed over all four of the detector chips. During the eclipse phase, all four chips 



were on but during the rest of the observation only one of the four chips was on, resulting in a 
~45% collection efficiency due to the placement of the source near the chip boundary. The spectra 
were extracted with the standard channel binning applied by XSELECT resulting in 512 energy 
channels. The spectra from the different chips of each detector and from the two SIS detectors 
from the same time interval were added using ADDASCASPEC. The energy channels in the range 
2.9-8.0 keV were binned additionally by a factor of two. Additional binning was applied to the 
energy channels in the range 8.0-10.0 keV so that each channel in the eclipse spectrum had at least 
50 counts. 

The analysis tools we used to compute the detector response include the reduction in effective 
area due to the fact the source areas we have chosen do not include all of the photons focussed 
by the mirrors. However, the background regions also contain some source photons which are 
subtracted from the source spectrum in computing the background selection region. Our source 
region (for all chips on) contains approximately 80% of the photons for a point source at the center 
and our background region contains approximately 15% (Serlemitsos et al. 1995). The background 
spectrum is multiplied by the area of the source region and divided by the area of the background 
region and then subtracted from the source spectrum. Because the background region is larger 
than the source region by a factor of 3, approximately 5% of the source flux is subtracted from the 
source spectrum. This effect is not accounted for by our analysis. This problem is compounded by 
the fact that the image of Cen X-3 is further smeared, in a way which depends on energy, due to an 
X-ray halo and this halo is delayed in time relative to the direct photons. However, the radius of 
the halo is approximately the size of the source region and the surface brightness due to the halo is 
generally no larger than that due to the image of the direct photons (Woo et al. 1994) and therefore 
leads to a similar subtraction of source photons. Because of the complexity of these effects, we do 
not try to account for them explicitly but note that fluxes (and quantities proportional to flux) 
which we derive are too small by approximately a factor 10%. This error does not qualitatively 
affect any of the results we derive in this work. 



3. Single-Zone Spectral Models 

Emission spectra from X-ray binaries are generally interpreted with the assumption that pho- 
toionization from the compact X-ray source is the dominant source of ionization in the circumstellar 
plasma. This is justified by the relative values of the luminosities of the X-ray sources and the den- 
sities and linear scales of the system. While it is straightforward to measure X-ray luminosities 
and, where orbital parameters are available from pulsar timing and optical spectroscopy, linear 
scales in X-ray binaries, determinations of the matter density that do not depend on the assump- 
tion of photoionization equilibrium are less reliable. Photoionization equilibrium has been inferred 
directly from observed X-ray recombination spectra (e.g., Liedahl & Paerels 1996). In the case of 
Cen X-3 however, there are no obvious spectral signatures of recombination dominance. The ob- 
served spectrum consists mainly of Lya transitions from hydrogen-like ions which, in principle, can 



be produced in plasmas where collisions dominate the ionization. We therefore attempt to fit the 
eclipse spectrum of Cen X-3 with spectral models of emission from collisionally ionized (coronal) 
plasma as well as photoionized plasmas. 

The X-ray spectrum observed from Cen X-3 includes at least the following components: direct 
emission from the neutron star, X-rays from the neutron star which have been scattered in the wind, 
continuum and line emission from the wind and other circumstellar material, and X-rays scattered 
from interstellar dust grains. All of these components contribute to the observed continuum, but 
only the wind and circumstellar material can emit lines. We found that, in general, it was possible 
to fit the continuum using two power laws, each absorbed by a different column density. The 
more highly absorbed power-law may correspond approximately to the X-rays from the neutron 
star viewed through some dense component of circumstellar material and the less absorbed power- 
law to the neutron star continuum scattered in the extended wind and by dust. However, since 
our primary goal is to extract and interpret the emission line spectrum, we do not attempt to 
constrain the parameters of these power laws in a manner that would force them to correspond to 
physical sources of continuum emission (c.f., Ebisawa et al. 1996). We interpret the lines, except 
for the 6.4 keV Fe Ka line as emission from the extended wind. These plasma emission models are 
described in detail below. The spectral model is, 

^(e) = e-(^)^Hi I jpj^(,) + /pi^^^^(e) + e-W^H2 [j^.^^5(, _ .^.^j + /pi2(e)]} (1) 

where the power laws. 



/ph(.) = K,u [j^) , (2) 

/plasma IS the plasma emission model, (T(e) is the interstellar absorption cross-section of Morrison 
&: McCammon (1983), N-m are the absorption column densities, /une is the line photon flux and 6 
is the Dirac delta function, used for the the 6.4 keV Fe Ka line complex. The energy of the Fe Ka 
line, eiine, was constrained to be in the range 6.3-6.5 keV. The Fe Ka line was the only fluorescent 
feature required to fit the data in any of our fits. Ebisawa et al. (1996) identified an emission 
feature at energy 1.25±0.04 keV with the 1.25 keV Mg I Ka fluorescent line. However, in all of our 
spectral models described here, this emission feature is fully accounted for by Ne X Ly/3 (1.21 keV). 

When the neutron star is eclipsed by the companion star, the X-rays from the extended wind 
only are observed. Therefore, for the purpose of testing basic plasma emission models, we fit only 
the spectrum from eclipse. For our spectral fits, we used the XSPEC spectral fitting program 
(vlO.O, Arnaud 1996), importing our own models for emission from photoionized plasmas. 



3.1. Collisional Ionization Equilibrium 

We fit the eclipse spectrum, using for /piasma the emission spectrum of an isothermal plasma 
in collisional ionization equilibrium (CIE). Also referred to as coronal equilibrium, this describes 
a situation where recombination is balanced by collisional ionization by electrons, and ionization 



by radiation is negligible. We use the MEKAL model (Mewe et al. 1996) contained in XSPEC 
for the CIE plasma emission. The three parameters of the MEKAL model are the temperature, 
metal abundance, and the normalization. The emission processes in CIE are determined by two- 
body interaction rates, so the normalization of the flux is proportional to the emission measure 
(E = J ri^dV) divided by the square of the distance to the source. The best fit parameters are 
shown in Table 1, and the spectrum is plotted in Figure 1. 

The MEKAL model, with best-fit parameters, which reproduces most of the the observed 
emission lines, also accounts for all of the continuum emission below ~4keV through bremsstrahlung 
radiation from the same gas. For parameters for which the MEKAL model accounts for the observed 
line emission, the model also accounts for all of the continuum emission below ~4 keV. However, 
coUisionally ionized gas must also scatter X-rays from the neutron star. In the ASCA band, optically 
thin electron scattering reproduces the continuum shape of the X-ray source with a fractional 
luminosity relative to the source spectrum equal to: 



L Att r 






dV, (3) 



where cjt is the Thomson cross-section and r is the distance from the compact source. We can 
estimate this fraction by taking the orbital separation a as the linear scale of the system and setting: 

4.77- 
r = a, n, = {E/V)^/^ V = —a\ (4) 

Then, 

^ « {UTry^/^a^E'/'a-'/' = 1.6 x lO'^ (5) 

If the neutron star has an intrinsic spectrum in the 2-10 keV band which is a power law with 
photon index 1.5, and luminosity ~ 10^^ ergs^^, then, using the system parameters of Table 2, the 
Compton scattered luminosity in the 2-10 keV band should be ~ 10^^ ergs~^. This corresponds to 
a power-law normalization of /pi = 1.6 x 10^^ s^^ cm^^ keV^^. However, the upper limit on /pn 
from the fit, 1.1 x 10~^s~^ cm^^keV^^, corresponds to a scattered luminosity of 7 x lO'^^ergs"^. 
In principle, it is possible that the Compton scattered continuum could be reduced relative to the 
thermal emission due to clumping of the wind. If the wind is clumped with volume filling factor /, 
then the density would have to be increased by a factor f~^''^ to preserve the emission measure. 
The new Compton scattered flux would then be changed by a factor /^'^ since its magnitude is 
proportional to nV. A very small filling factor (/ ~ 10~^) would be necessary however and we 
regard this as unlikely. We therefore reject this model. 

Increasing the metal abundance in the gas increases the flux of the lines relative to the 
bremsstrahlung continuum, thereby providing "room" for a scattered continuum component. There- 
fore, we tried another fit in which we allowed the metal abundance to be free. For this fit, we tied 
the normalization of the first power law (pll) to the emission measure of the MEKAL component 
as described above (Eq. 5). XSPEC does not allow the user to set one parameter to the square 
root of another parameter, so we fixed the power-law normalization according to a trial value of the 



emission measure, fit, and iterated to find a best fit emission measure with a consistent power-law 
normalization. This procedure resulted in a statistically acceptable fit (Table 1, Figure 2). The 
lower limit on the abundance derived from this procedure, however, is 25 times solar, which we 
believe to be implausibly high, and so we reject this model as well. 



3.2. Photoionization Equilibrium 

Having rejected plausible collisional ionization spectral models, we proceed to fit the observed 
spectra using model emission spectra from photoionized plasmas. The radiation spectrum due to 
recombination consists of radiative recombination continua (RRC) from free electrons recombining 
directly to bound states and lines which are produced through radiative cascades. The volume 
emissivity of the recombination feature k is jk = ricfii^iakiT) where nj+i is the density of the 
recombining ion to which the electron recombines. The function a^iT) is an effective recombination 
rate coefficient (the recombination coefficient for recombinations that produce transition k) that 
depends on atomic parameters and the level population kinetics. In a gas which is optically thin 
and has density low enough such that all ions may be assumed to be in the ground state (i.e., 
the density is low compared to the critical density for all important transitions), the ion fractions 
and temperature are, for a given spectrum of ionizing radiation, functions only of the ionization 
parameter (^ = L/nr'^, Tarter, Tucker, Sz Salpeter 1969). Therefore, for a photoionized plasma, 
the emissivity of any feature may be written j^ = nlfk{S,), and, for the entire spectrum, ju/n^ is a 
function only of ^. 

To calculate the temperature and ion fractions as a function of ^, we need to know the X-ray 
spectrum of the neutron star. The best way to do this is to observe the X-ray spectrum from 
the neutron star directly while outside of eclipse. This ASCA data set includes data outside of 
eclipse. However, the observed spectra in these phases outside of eclipse differs greatly from the 
intrinsic spectrum of the neutron star because at the time of the observation the line of sight to 
the neutron star contained circumstellar material that substantially absorbed the direct X-rays. 
When circumstellar material passes in front of the neutron star, the direct X-rays are absorbed 
preferentially at low energies. Even when direct low energy X-rays are absorbed completely, low 
energy X-rays are still observed from scattering and emission from the extended wind, just as during 
eclipse by the companion star. The residual X-rays from scattering are generally not pulsed. The 
X-rays from Cen X-3 in this data set show both a spectrum which is deficient in low-energy X-rays 
compared to other observations (Nagase et al. 1992; Santangelo et al. 1998) and also are not pulsed 
at low energy (Ebisawa et al. 1996). Nagase et al. (1992) used similar evidence to demonstrate 
that a "pre-eclipse dip" observed by Ginga was due to absorption. In these out-of-eclipse intervals, 
the large absorption of the direct X-rays makes the strength of the direct fiux comparable to and 
difficult to distinguish from the residual flux from scattering and wind emission, and it is impossible 
to isolate the flux from the neutron star in the observed spectrum. For the intrinsic spectrum of 
the neutron star, we instead use a spectrum obtained by Santangelo et al. (1998) with BeppoSAX. 



Burderi et al. (2000) have shown that in this observation, the X-ray flux is pulsed down to the energy 
band below 1.8 keV indicating that it originates in the proximity of the neutron star. There is an 
additional advantage to using BeppoSAX in that it can measure the spectrum up to ~ 200 keV. For 
C ^ 10^1 the ions are nearly fully stripped, and the temperature is determined by Comptonization. 
The Compton temperature is highly sensitive to the presence of high energy photons. We used the 
"Lorenzian" model of Santangelo et al. (1998), extending the power-law continuum to low energy 
by setting the absorption column to zero, and excluding the cyclotron absorption feature, since a 
value for the depth of this feature is not provided. 

We used the XSTAR program (vl.43, Kallman & Krolik 1999) to compute ion fractions and 
temperatures for 100 uniformly spaced values of log^ from —2 to 6 using the Santangelo et al. 
(1998) spectrum. In the XSTAR calculation, we used a luminosity of lO^'^ergs"^ and a constant 
density 10~^ cm~^ — again, for an optically thin gas, the absolute values of the luminosity and the 
density are not important. We then calculated the recombination spectrum at each value of ^ using 
the ion fractions and temperatures from XSTAR, and line and RRC powers calculated using the 
Hebrew University/Lawrence Livermore Atomic Code (HULLAC, Klapisch et al. 1977) and the 
Photoionization Cross-Section code (PIC, Saloman, Hubble, & Scofield 1988) for recombination 
cross-sections. These atomic models have been used to analyze the recombination spectrum of 
Cyg X-3 (Liedahl & Paerels 1996) and Vela X-1 (Sako et al. 1999) and are described further there. 

We fit the eclipse spectrum using for /piasma our recombination model spectra plus a cor- 
responding bremsstrahlung component. The bremsstrahlung component is the BREMSS model 
from XSPEC. We tied the normalization of the bremsstrahlung component to that of the recom- 
bination component such that the emission measures would be equal. We set the temperature of 
the bremsstrahlung model to be a locally linear approximation to the function T(log^) computed 
with XSTAR in the neighborhood of the best fit value of log^. The best fit spectrum is shown in 
Figure 3, and the best fit parameters are shown in Table 1. The best fit value of the ionization 
parameter (log^ = 3.19) corresponds to an electron temperature of 0.41 keV. At this temperature, 
the radiative recombination continua are broad, and so it is not surprising that they do not appear 
prominently in the spectrum. It is because of this that Ebisawa et al. (1996) were able to fit the 
spectrum using only lines and not RRC. This demonstrates that though the presence of prominent 
narrow RRC in emission spectra indicates photoionization (Liedahl & Paerels 1996), the lack of 
narrow RRC does not necessarily indicate some other type of equilibrium. It can be seen that the 
bremsstrahlung component is relatively faint, thus allowing a power-law continuum at low energies. 
The luminosity implied for the power-law component #1 is 4 x 10^^ ergs^^which is much closer to 
the expected Compton scattered luminosity of rsj lO^^ergs"^. A quasi-continuum of blended lines 
and RRC accounts for 21% of the photon flux in the 1-3 keV band. This differs significantly from 
the individual line fits of Ebisawa et al. (1996) in which the lines constitute only 5% of the flux in 
this band. 

From this single-zone fit, it is possible to make a crude estimate of the X-ray luminosity and 
parameters of the wind in Cen X-3, if we take the orbital separation to be the characteristic linear 



10 



scale of the system, and use our best fit value of the emission measure as in Equation 4. We then 
we find for a characteristic wind density n --^ 5 x 10^*^ cm^"^((i/10kpc)~^. From the definition of ^, 
we then have 

L = S,nr'^ ~ 1.4 x 10^^ ergs~^((i/10kpc)~^ (6) 

Note that this value of the luminosity is derived only from the flux of the recombination spectral 
features and the known system dimensions, yet it comes very close to the luminosity derived for 
Cen X-3 by measuring the continuum directly during the high state (Nagase et al. 1992; Santangelo 
et al. 1998). Ebisawa et al. (1996) derived similar, though somewhat larger, estimates for the 
characteristic density and luminosity (1.6 x 10^^ cm^^ and 1.8 x 10^^ ergs"^ for (i=10kpc) using an 
emission measure derived from the fluxes of the 6.7 keV and 6.9 keV lines of, respectively, helium- 
like and hydrogen-like iron in the eclipse and egress spectra and by assuming a value of log .^ = 3.4 
from their ratio. Our values differ because we use slightly different values for the linear scale and 
also because our values of S, and E are derived from fitting all of the lines and not just the hydrogen- 
and helium-like iron lines. 

We now estimate the characteristics of the stellar wind from this derived value of the density. 
In a spherically symmetric steady-state wind, continuity demands that: 

^ ATrR^fin{R) (7) 



v{R) 



where M is the mass-loss rate, v is the velocity, R is the distance from the center of the star, 
n is density of hydrogen (neutral and ionized), and the quantity fi is the gas mass per hy- 
drogen atom which we take to be lAmn. Again using the fiducial dimensions, M/v ~ 4 x 
lO^^M0yr~^ (kuis^^)^^{d/10kpc)^^ which, as will be discussed in § 6 below, is near values for 
isolated stars with spectral type similar to V779 Cen. While this model provides a statistically 
acceptable fit to the data, it is unlikely that the entire emission line region can be characterized by 
a single value of ^. Therefore, we suggest that this model provides no more than a semi-quantitative 
description of the data. In the next section we use explicit wind models and spectra of all of the 
full phases of the observation to derive these wind parameters more accurately and to explore the 
wind geometry. 



4. Global Photoionized Wind Models 

We use the recombination spectra described in the previous section to derive spectra for explicit 
wind density distributions. In order to calculate the global recombination spectrum of photoionized 
gas, we use the differential emission measure formalism (DEM, Sako et al. 1999). For an arbitrary 
distribution of gas, the total spectrum of recombination radiation is given by 



, _ , MO 

J-'v 



dE 
d\ogi 



dloge, (8) 
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where the quantity in brackets is the DEM distribution, which hereafter will be referred to as D{(^). 
In practice, we integrate only the volume of gas that is visible to the observer and refer to this as 
the apparent DEM. To calculate the DEM distribution, it is necessary to know only the density 
distribution of the gas and the luminosity of the radiation source. To calculate ji/(^)/ng, it is 
necessary to know only the spectral shape of the ionizing radiation and the elemental composition 
of the gas. In this work, we calculate the spectrum of recombination emission from the wind of 
Cen X-3 for different models of the wind using different luminosities for the neutron star, but 
always using the same spectral shape for the neutron star X-ray emission and always assuming 
solar abundances (Anders & Grevesse 1989). The DEM formalism requires us to evaluate JuiO/^e 
only once for a number of values ^, and then, for different matter distributions and luminosities, 
to calculate the DEM distribution and evaluate the integral in Equation 8. 

To calculate the DEM distribution for a given density distribution and luminosity, we divide 
the binary system into spatial cells. We then calculate ^ and the emission measure for each spatial 
cell. The emission measure for each cell is then added to a running sum of the emission measure 
for a ^ bin. Dividing the emission measure in each bin by the width of that bin in log ^ gives the 
differential emission measure. If the density distribution is described by parameters such that the 
density everywhere scales linearly with one parameter, which we will call rj, then to calculate the 
DEM distribution for new values of ij and the luminosity L it is not necessary to recalculate the 
DEM distribution by summation over the spatial cells. For new values L' and r]' , the new DEM 
distribution D' is related to the old DEM distribution D by 

This identity is derived in Appendix A. For a change in the luminosity and density parameter such 
that L' /rj' = L/t], the DEM distribution is changed only by a constant factor (rj' /i])'^ = {L' /L)'^ 
and therefore so is the total emission spectrum. 



4.1. DEM Distributions for Model Winds 

As alluded to in § 1, the structure of stellar winds in X-ray binaries are likely to be rather 
complicated. In addition to the effects of X-ray photoionization on the radiation driving, the 
gravity of the compact object and orbital motion may have a significant effect on the wind density 
distribution (e.g.. Friend & Castor 1982). Though numerical simulations may account for many of 
these effects, we choose instead to use a simple, spherically symmetric density distribution with free 
parameters such that the density is easily recalculated for a change in parameters. We approximated 
the geometry of the X-ray binary system as a spherical star in orbit with a point-like neutron star 
emitting X-rays isotropically. For the dimensions of the system, we assumed the values in Table 2. 
We used density distributions for spherically symmetric, radiation driven winds (Castor et al. 1975; 
Kudritzki et al. 1989) from the companion star. 
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Modeling accretion disk winds introduces a new set of complications. They may be driven 
thermally (e.g., Woods et al. 1996), radiatively (e.g., Proga, Stone, & Drew 1999), or flung out 
along rotating magnetic field lines (magnetohydrodynamically, e.g., Blandford & Payne 1982; Proga 
2000). While stellar winds may be symmetrical in two dimensions (azimuth and altitude), disk 
winds are necessarily symmetric in no more than one dimension (azimuth). However, in order 
to derive some observational constraints on a disk wind hypothesis, we use the same spherically 
symmetric symmetric radiation-driven velocity profile as for the stellar wind but center the wind 
on the neutron star. 

The explicit form of the wind velocity profile is given by: 

viR)=vo + iv^-vo)il-Rin/Rf (10) 

where vq and foo are the wind velocities at the stellar surface and at infinity and the parameter 
P describes the acceleration of the wind. With /? = 0, this equation describes a wind which is 
immediately accelerated to its terminal velocity. The variable R is the distance from the center of 
the wind (the center of the companion for the stellar wind or the neutron star for the disk wind) . 
For the stellar wind, the wind begins at the surface of the companion so i?in = -R*. For a disk wind 
however, no such natural inner radius exists and therefore, for the disk wind, we make i?in an extra 
parameter of the model and assume that the volume for which R < R[^ is empty. The value of i?in 
may correspond, for example, to a characteristic radius on the accretion disk from which the wind 
arises. Because an accretion disk may be no larger than the Roche lobe of the accreting object, the 
Roche lobe radius should provide an approximate upper limit on R^^. 

A rearrangement of Equation 7 gives the density distribution. 

In order to make use of the scaling relation of Equation 9, we include the velocity profile explicitly 
(Equation 10) and reexpress the density as follows: 

n{R)= I — ] {A-K^i)-^R-^ 



V, 



oo 



-+fi--Ui-— '' 

t'oo V Woo/ V R 



(12) 



According to this parameterization, the parameter M/voo plays the role of rj. The shape of the 
DEM distribution is then then a function of the parameters Lvoo/M, /3, vq/voo, and, for the disk 
wind, -Rin- For a given choice of those parameters, the magnitude of the DEM distribution is 
proportional to (M/foo)^ or, equivalently, L^. 



4.2. Calculation of DEM Distributions and Global Spectra 

The wind distributions we use are symmetric to rotation around the line containing the neu- 
tron star and the center of the companion star. This allowed us to conserve computational re- 
sources by using rings around this axis as the spatial cells to calculate the DEM distributions. 
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The exclusion of the region occulted by the companion star from the apparent DEM breaks this 
rotational symmetry but we were still able to conserve computation, using rings as the spatial 
cells, by computing the fraction of each ring not occulted by the companion star and multi- 
plying the emission measure for each ring by that factor. We used the binary system param- 
eters from Table 2, and fixed vq/voo at 0.015 which corresponds approximately to the photo- 
spheric thermal velocity of an 0-type star for vq if Wqo is of order 1000 km s^^. For the spec- 
tra which were accumulated during an orbital phase interval longer than 0.02, we averaged the 
apparent DEM over phases separated by no more than that. Our parameter grid contains the 
values (3 = (0.0,0.1,0.2, ..., 1.5), \og{Lvoo/M) = (-0.5, -0.33, -0.17,0,0.17, ..., 2.5) in units of 
-|^Q37gpgg-i £qj, j^^ 1000 km s^^ for v^o and IQ^^ MqJI^'^ for M. For the disk wind, the grid also 
contained R\a = (3.16, 5.62, 10, 17.8, 31.6) Rq. The resultant spectra were stored in a FITS format 
XSPEC ATABLE file (Arnaud 1995). 

In order to demonstrate the appearance and behavior of the DEM distributions, we plot DEM 
distributions and contour maps of log ^ using the best fit parameters found in the following section 
(§ 4.3, Table 3). In Figure 4 we show a map of log^ for the companion star wind and in Figure 5 
we plot the apparent DEM distributions for orbital phases representative of the observed phases. 
In Figures 6 & 7 we show the same plots for the disk wind model. However, as mentioned in § 4.1, 
we expect the Roche lobe radius to be an upper limit on R\^ and so we use R\^ = 3.4i?0, the Roche 
lobe radius as determined from the parameters of Table 2 and the formula of Eggleton (1983), for 
these plots instead of a best fit value of i?in. In Figure 7, it can be seen that the apparent emission 
measure is dramatically reduced during eclipse for a disk wind inner radius significantly smaller 
than the companion star radius. 



4.3. Spectral Fitting 

For both the disk and the stellar wind models, we simultaneously fit the spectra at all four 
phases, using for /piasma the recombination spectrum for the observable wind at each phase as 
discussed above. We allowed the parameters of all of the spectral components except for /piasma 
to vary independently for the four observed phases. We show the best fit values and associated 
errors in Table 3. We show the spectral fits for the stellar wind in Figure 8. The best fit model 
spectrum for the eclipse phase is shown at high resolution in Figure 9. Since our focus is on the line 
spectrum, we let the continuum parameters vary to fit the continuum shape without constraining 
them to physically meaningful bounds. Some strange results were obtained, such as the very large 
columns, normalizations, and photon indexes for the second power law in the egress phase (for 
both the stellar and disk wind models). Though the normalizations are very large, the fluxes of 
this component are comparable to that in the other phases. We note that our best fit value of 
the neutron star luminosity from the stellar wind model is very near the value obtained by direct 
measurement of the unocculted continuum. We emphasize that this luminosity value is obtained 
only from our recombination spectra for wind models, and is determined almost entirely by the 
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observed line fluxes, not from any direct measurement of the broad band flux. 

Both wind models provide acceptable values of x^ ) though for the disk wind the fits favor values 
of the inner radius which are at least comparable to the radius of the companion star. In Figure 10, 
we plot the value of x^ for the fits as a function of the assumed inner radius. As mentioned in 
§ 4.2, values of Ri^ much smaller than the companion star radius lead to dramatic variations in 
the DEM distribution, and therefore the emission line fiux over eclipse. The observed emission line 
spectrum, however, varies only moderately across the eclipse. As the preferred radii are very large 
compared to the neutron star Roche-lobe radius, we consider the possibility that the wind could 
be dominated by matter arising from the accretion disk unlikely. 

We compare the lines fiuxes predicted by our best fit stellar wind model with line fluxes 
measured by Ebisawa et al. (1996) in Figure 11. According to our model, the lines from elements 
with lower values of Z vary somewhat more as a function of orbital phase. This is because the lines 
from low-Z ions are produced more efficiently in the region near the surface of the companion star, 
where the ionization parameter is lower, and which is occulted more during eclipse. Unfortunately, 
the quality of the current data is not sufficient to test this prediction. Ebisawa et al. (1996) noted 
that the H-like and He-like lines of iron were undetectable during the first phase ( "pre-eclipse" ) 
from which they inferred that the overall ionization parameter was lower at that phase. Our 
analysis presumes that we observe the same symmetric wind, only occulted differently, at every 
phase. Therefore, in our model, the line fluxes must be the same at symmetric phases. However, 
our predicted line fluxes are greater only by a factor of a few from their la limits and so our 
assumption that the wind density distribution and luminosity are constant is not excluded. 



5. Optical Depth of the Stellar Wind 

Until now, we have assumed that the stellar wind is optically thin to X-rays from the neutron 
star. For the densities and luminosities of our best fit model, the state of X-ray photoionized 
gas is affected little by optical depth inside an ionized Stromgren-type zone, where the charge 
state distribution of helium is dominated by He^'^ but outside this zone, where the charge state 
distribution of helium becomes dominated by He"*", the effects of optical depth quickly become 
important (e.g., Kallman & McCray 1982). The size of a Stromgren zone is usually derived for 
spherically symmetric nebulae (e.g., Osterbrock 1989) with the condition that the rate at which the 
central source emits photons capable of ionizing helium is equal to the total rate of recombinations 
which do not produce ionizing photons (i.e., recombinations to excited states). However, the 
Stromgren radius may also be estimated in a nebula which is not spherically symmetric for any 
direction from the ionizing source using the condition that the rate of helium ionizing photons 
emitted per solid angle is equal to the total rate of recombinations per solid angle: 

1 f°° L f^^ 

—— I T-dl' = / aB^e^He++'^e^ dr (13) 
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where ob is the recombination rate to excited states and uq is ionization threshold of He^"*" (hi^o = 
54.4 eV). If we assume that the spectrum of Cen X-3 is a power-law with photon index 1, cut off 
at energies above E'cut = 15 keV (approximately the Santangelo et al. 1998 spectrum) then the 
integral on the left hand side is 

LE-^, ln(i?eut//i^o) = ^;^^ = 2-3 x 10^^ s'^ (14) 

for L = 10^^ ergs~^. The computation of the integral on the right-hand side is complicated by the 
fact that the recombination coefficient, ob, is a function of temperature. However, the temperature 
dependence of ob is only --^ T^^'^ and so we can get a good estimate by assuming that the 
nebula is isothermal. In the calculations of Kallman &; McCray (1982), the temperature just 
inside of the He'''^/He"'" boundary is approximately 10^ K and so we use for ob the constant value 
2 X 10~^^ cm^ s^^ (Osterbrock 1989, Table 2.8). In the ionized zone, nearly all of the helium is 
in the form He"'""'' and so nHc++ = 8.3 x lO'^nc almost exactly. The condition that the nebula is 
optically thin becomes: 

/ n^r^dr < 1.1 X 10^^ cm"^ (15) 

Jo 
We note that the quantity on the left hand side is the emission measure per solid angle. For our 
wind, from the neutron star to the face of the companion along the line of centers, this integral is 
2.72 X 10^'^ cm~^ and in the direction away from the companion to infinity, it is 2.98 x 10^^ cm~'^. 
Therefore, it is clear that the assumption that our model wind is optically thin, which we have used 
to calculate emission spectra, is self-consistent. 



6. Discussion 

We note that our results for the global wind models can be generalized beyond the homoge- 
neous, spherically symmetric distributions we have used. If the wind is clumped on a scale which is 
small compared to the system such that the ionization parameter inside the clumps and the total 
emission measure is the same as in a corresponding homogeneous model, then the observed spec- 
trum will be unchanged. To construct such a model where the clumps have a filling factor /, the 
density must be increased by a factor /^^' ^ to keep f n^dV unchanged and so to keep ^ unchanged, 
L must also be increased by a factor /~^'^. Because our best fit luminosity is already at the high 
end of the range of luminosities observed for Cen X-3, the wind cannot be clumped with a filling 
factor significantly smaller than unity. A similar, though no longer exact, extension can be made 
to winds that, instead of being confined in small clumps, are confined in solid angle. An example 
of this scenario is a wind which exists on the surface a cone. Such a geometry is believed to exist 
in broad absorption line quasars and may be due to preferential launching of the wind along rays 
from the central object which graze the surface of the accretion disk (e.g., Murray et al. 1995; 
Proga 2000). For such a geometry, the factor VL/Ati, where VL is the solid angle to which the wind 
is confined, plays the same role as / above. For the same reasons then, we can exclude these types 
of winds. 
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Both the stellar wind and the disk wind models give statistically satisfactory fits to the data. 
However, the best fit values of the inner radius in the disk model are very large compared to the size 
of the accretion disk. If the extended reprocessing material originates primarily from an accretion 
disk wind, a geometry in which the wind is prevented from producing recombination radiation inside 
a radius which is approximately equal to companion star radius (which is ~4 times the maximum 
size of the accretion disk) is required. We believe that it would require fine-tuned conditions to 
produce such a geometry, and we consider the possibility that a disk wind dominates the extended 
circumstellar material unlikely. 

As we have noted, in our analysis, the mass loss rate and terminal velocity can be obtained 
only in the combination M /voo- If Doppler shifts of lines due to motion of the wind could be 
measured, these two parameters could be determined independently. The errors in the line energies 
quoted by Ebisawa et al. (1996) correspond to velocity upper limits of ~2000kms~^ for the iron 
lines and ~6000kms~^ for the lower energy lines. Therefore, only an upper limit on the mass loss 
rate of the highly ionized wind of ~ 3 x 10~^ Mq yr~^ can be obtained. 

The fact that the value of the luminosity that we obtain from our spectral line fits is very 
close to the value determined from direct measurements of the continuum gives us confidence that 
the parameters of our model correspond to unique physical values that accurately describe the 
stellar wind. For the stellar wind, the best fit values of M /voo and (3 are similar to those found in 
isolated stars. For example, for the 06.5 laf star HD 163758 and the 06.5 Ill(f) star HD 190864, the 
values of M are 6 x 10~^ and 1.5 x 10~^ Mqyt~^ and the values of Voo are 2200 and 2500 km s~^ 
(Lamers et al. 1999) which corresponds to values of Mjvoo of 2.7x10"^ M0yr~^ (kms^^)^-*^ and 
0.6xl0"9^^yj.-i(knis"i)"i as compared with 1.56±0.12 x 10"^ M© yr^^ (kms^^)"! (Table 3) 
for our fits to the Cen X-3 spectrum. The values of /3 for these two isolated O stars are 0.7 and 
0.8 (Lamers et al. 1999), which are also near our derived value of 0.57lgQ7 (Table 3). Therefore, in 
our model, the best fit wind parameters are roughly consistent with those of the normal radiatively 
driven winds in isolated massive stars. Considering that the wind is highly ionized and the radiation 
driving mechanism which governs the structure of the wind in isolated massive stars cannot function, 
this result is very surprising. 

To within an order of magnitude or so, the total emission measure in a smooth wind is given by 
E ~ [M / Voo)^ / {fniR^,) . Therefore, to the extent that the companion stars in HMXBs are roughly 
the same size, lose mass at roughly the same rate, and produce smooth winds with roughly the 
same terminal velocities, the total emission measures should be approximately equal. Numerically, 
the above estimate gives, for typical parameters, E ~ 10^^ cm~^. (Note that this estimate applies 
equally well to isolated early-type stars.) Clumping in the wind, however, will tend to increase 
the emission measure. On the other hand, clumping also tends to decrease the local value of ^, in 
effect, removing gas from the high-^ end of the DEM distribution and adding gas to the lower-^ end. 
Thus clumping reduces the wind's X-ray recombination line emission, and increases its fluorescence 
emission. 
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We have argued that the Cen X-3 wind is smooth. Sako et al. (1999) have shown that the 
wind in Vela X-1 is highly clumped. In Figure 12, we compare the DEM distribution derived here 
for Cen X-3 with that derived, also using X-ray recombination, by Sako et al. (1999) for Vela X-1. 
The value of M /voo we have derived here for the wind of Cen X-3 is comparable to the value for 
Vela X-1 derived by a number of methods (Dupree et al. 1980; Sato et al. 1986; Sako et al. 1999). As 
one caveat, we point out that for Vela X-1 the contribution from fluorescing material, which would 
cause an upturn in its DEM curve below log^ ~ 1.5, is not included in the plot.^ Above log.^ = 2.0, 
however, the differences are real, and are quite striking. The relatively small DEM magnitudes in 
Vela X-1 are a consequence of the fact that most of the mass in the wind is "locked up" in clumps 
of high density, hence, low ^. For example, near log^ = 3.0, where lines from He-like and H-like 
iron are produced, the DEM magnitudes vary by ~ 100, which means that these lines are ~ 100 
times more luminous in Cen X-3 than in Vela X-1. This does not translate into excessively large 
line equivalent widths in Cen X-3, however, since its X-ray luminosity is also roughly 100 times 
higher. 

With this comparison arises the question as to what conditions are required in order to produce 
and/or destroy clumps in X-ray irradiated winds, a subject beyond the scope of this paper, but 
possibly one that bears on the nature of the wind driving force. The absence of substantial clumping 
in Cen X-3, as inferred from the modest phase variations of the iron Ka equivalent width, is 
demonstrated by, and is consistent with, the large DEM magnitudes for log^ > 2. It has been 
suggested (Sako et al. 1999) that the existence of a clumped wind component in Vela X-1 allows 
normal radiative driving by the UV field of the companion. Clearly, that is disallowed in the case 
of Cen X-3. We therefore conclude that the wind is most likely driven by X-ray heating of the 
illuminated surface of the companion star as proposed by Day & Stevens (1993). 

We look forward to Chandra observations of this system. High resolution spectroscopic data 
will allow us to measure independently the lines of the helium-like 2^1 triplets and therefore 
determine the ionization mechanism directly (Liedahl 1999). We will also be able to measure much 
smaller velocities (of order 100 kms^^) and therefore be able to set much tighter constraints on the 
mass-loss rate of the companion and the dynamics of the wind. 



A. Differential Emission Measure 

We have defined the emission measure as fri^dV. We now define the ionization parameter 
limited emission measure as: 

m) ^ / nldV, (Al) 



^Mapping the DEM distribution for low-^ material based upon X-ray spectra requires much higher spectral 
resolution, so that fluorescent line complexes can be resolved into their respective charge states. 
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i.e., the emission measure in the volume in which the ionization parameter is less than some value. 
The linear differential emission measure is then 



This infinitesimal volume is the space between two surfaces, each of which is defined by a single value 
of the ionization parameter. The infinitesimal distance between the surfaces is (5^/|V^|. Therefore, 

'"^^^-f "' 4S^ (A3) 



di Js(() ivei 

where S{S,) and dS specify the integral over the surface specified by the ionization parameter ^. 
We define the differential emission measure 

DiO^^=[ T^dS. (A4) 

^^^ dlogC Js^^) iVlog^l ^ ^ 

Suppose we have a model for the density distribution such that the density scales linearly with 
the parameter r]. Consider a change of parameters ij ^ rj' and L ^ L' . The surface defined by ^ 
for the unprimed parameters is the same surface defined by 

«'-0e (A5) 

Therefore, D'{S,') can be related to D{^) by transforming the integrand in Equation A4. Under 
the parameter change, log^ differs from log^' by an additive constant and so |Vlog^| = |Vlog^'|. 
However, n^ is modified by the factor (f]' /i])'^- Therefore, 



e 

2 



D'ia = ( ^ ) DiC). (A6) 



This is Equation 9. 
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Table 1. Fit Parameters for Eclipse Spectrum 





Spectral Line Model 






Variable 








Abundance 




Parameter 


MEKAL 


MEKAL 


Photo 


A^Hl (1022(,j^-2) 


n qfi+o-03 


0.91±0.02 


1.03±0.04 


E {W^^cm^^){d/10 kpc)-2 


o /^o+0.06 
•-••"^-0.08 


Q .7+0.09 


2.58^°:- 


logC 






3.19tH? 


kT (keV)^ 


7 ^+0-4 


8.9±0.5 


0.41±0.01 


Z/Zq 


1 (fixed) 


391^- 


1 (fixed) 


ai 


1.5(fixed) 


1.5(fixed) 


1.50±0.09 


/pli (10-3 s-^cm-^keV-i) 


<0.11 


5.6(fixed) 


6.7±0.6 


nm (1022 cm-2) 


71+10 


74±14 


^otll 


eiine (keV) 


a ^1+0.018 

'3-'±-L_o.o20 


6.393+Hi? 


6.388^°:°?° 


/iine(10-^phcm-2 g-i) 


q O+0.9 
-3-0-0.8 


A 8+1-4 


q q + 1.2 


02 


1.5±0.4 


1 fi+0-5 
J^-O-0.6 


1 5+°-'^ 

-L-^-0.4 


Ipl2 (10-3phs-icm-2) 


-^0+57 
'3y-25 


11+1^7 


-,7+79 


xVd.o.L 


176/170 


146/170 


176/169 


probability 


35% 


90% 


33% 



^In the case of the photoionized model, kT is an explicit function of 
log^. The error on the temperature is estimated using a local linear 
approximation to T(log^). 

Null hypothesis probability 
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Table 2. Assumed System Parameters for Cen X-3/V779 Cen 

Parameter Value Reference 

i?* (companion radius) 11.8 Rq a 

a(orbital separation) 19.2 Rq a 

i(inclination angle) 70.2° a 

Distance 10 kpc b 



References. — (a) Ash et al. 1999 (b) Hutch- 
ings et al. 1979 



Table 3. Global Wind Model Fit Parameters 







Stellar 


Wind 






Disk Wind 






Parameter 


Pre-cclipse 


Ingress 


Eclipse 


Egress 


Pre-eclipse 


Ingress 


Eclipse 


Egress 




Nni (1022cm-2) 


0.95±0.08 


„ 07+0.02 


1 nn+0-02 

^■UU_o.04 


o.94«;SI 


0.86±0.08 


n 04+0.05 


1.03±0.04 


n 07+0.07 




ai 


o.59t°:lI 


n Qfi+0.09 


1 49+0.11 


-0.181HI 




n Qfi+010 


1 -1+0.12 
1-41-0.05 


-0.191H? 




Ipll (10-3phcm-2 s-ikcV-i) 


17.4t^;« 


7 q+0.8 

'■•^-0.5 


5.8t^-i 


3.lt^i 


179+2.5 


7.3t«:^ 


5 5+0.7 

0.0_g 2 


3.09l«i? 




Nii2 (10^2 cm-2) 


35l| 


140tl^ 


50ti! 


183li^ 


38t? 


154tf6 


46t\l 


197+?i 




eiine (keV) 


fi 074+0.023 


6.3ntO:S« 


6.39linil 


6.390±0.020 


R ovR+0.024 
D.ci/-D_g Q28 


f. 071+0.017 

O-^' ^--0.022 


6.392±ni5 


6.39ltn^? 




/line(10-*phcm-2 s-1) 


23±5 


40if3 


q+0.7 

•^■•^-0.6 


2.2+J-^xl02 


24±5 


49±1« 


9+0.8 


271^2 




«2 


i-5tS:l 


q+0.6 

■^■•^-0.9 


1 ^+0* 

^■^-0.5 


7 5+0.9 


1 fi+0-5 

^-■"-0.4 


3.6±?:g 


1-51EJ1 


8.0+?i 




-fpl2 (10-^phcm-2s-lkcV-i) 


2.7±«xl02 


6.1±*«4X103 


20tll 


8.3^620 yxio^ 


3.3t5;|xl02 


lAtflxm 


3 10+39 

18-13 


2.6^190^x108 




Rm(i?0) 














>15 






/3 




0.571H? 






n 51+0.17 

U-Ol_()07 






Lv^/M "^ 




6.8l?:° 








15^22 






M/i,^[d/(10kpc)]-i- 




1.56±0.12 






1 3+°-2 
J^-^-0.6 




1 


L[d/(10 kpc)]-i ^ 




10.7 


+2.9 
-2.4 








20±9 




K3 

CJl 


X^ 




552/684 






550/683'' 




1 



''These three expressions are combinations of only two non-degenerate model parameters. The units for these values arc such that M has units of 10 Mq yr 
units of lOOOkms"^ and L units of lO^'^ ergs^^. 

''for i?in=30 Rq 
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Table 4. 


Stellar Wind Model 


n = 


2 -^ 


1 Line Fluxes'^ 






Line 

Energy(keV)^ 






Phase Range 




Ion 


-0.31 --0.29 


-0.23 - -1 


0.08 


-0.08 - 0.13 


0.14-0.20 


NeX 


1.022 


4.7 




2.7 




1.9 


3.0 


MgXII 


1.472 


1.7 




1.0 




0.7 


1.0 


Si XIV 


2.006 


2.0 




1.3 




1.0 


1.4 


SXVI 


2.621 


1.7 




1.2 




0.9 


1.2 


FeXXV 


6.667 


4.9 




3.6 




3.2 


3.6 


Fe XXVI 


6.966 


2.8 




2.2 




2.0 


2.3 



^In units of 10 ^ ph cm ^ s ^ are summed over line complexes and energies are for 
centroids of complexes. 

'^For He-like ions, the centroid energy of the line complex. 
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Fig. 1. — The fit to the eclipse spectrum using the MEKAL cohisionally ionized plasma model. 
The crosses represent the data. The solid line represents the total model spectrum and the dashed 
and dotted lines represent the components of the model spectrum. The MEKAL model (lines and 
bremsstrahlung continuum) accounts for nearly all of the observed flux below ~5 keV. Power law 
2 and the 6.4 keV line make up the remainer of the flux at high energy. No room for a Compton 
scattered continuum remains below 5keV. 

Fig. 2. — The eclipse spectrum fit with the MEKAL model with variable abundance. Here, the 
power law is larger than the MEKAL component at low energies though the emission features are 
approximately the same as in solar abundance fit. 

Fig. 3. — The fit using the single-zone recombination model. The recombination line and RRC 
spectrum (second largest component in 1-2 keV region) is plotted as a component separate from 
the 0.4 keV bremsstrahlung (small component peaking near 1.3 keV). In PIE, the observed lines are 
produced at a much cooler temperature than in CIE. Therefore, in PIE the bremsstrahlung contin- 
uum is much weaker and there is room for a Compton scattered continuum (which is represented 
by the two large continuum components). 

Fig. 4. — Contour plot of ionization parameter distribution for the stellar wind in the orbital plane. 
Contours of constant log.^ are indicated. Dashed lines delineate the region which is not visible when 
the system is viewed from an orientation opposite the neutron star. The position of the neutron 
star is indicated by crosshairs. 

Fig. 5. — Apparent DEM distributions for the best fit stellar wind model parameters as viewed 
at three phases characteristic of the ASCA observation. At eclipse center (phase 0) the regions of 
highest ^ (near the neutron star) and the regions of lowest ionization parameter (near the stellar 
surface) are highly occulted. Away from eclipse center the DEM increases at all values of ^ as more 
of the wind becomes visible, though the increase is most dramatic at the highest and lowest values 
of ^. 

Fig. 6. — Contour plot of ionization parameter distribution for a disk wind with parameters of 
Table 3 except with R^^ = SARq, the Roche-lobe radius. In this model, the region inside the 
innermost contour has zero density. At large radii, the density varies as R~^ and so ^ approaches 
a constant value. 

Fig. 7. — The apparent differential emission measure for the disk wind for orbital phases character- 
istic of this observation. The model parameters are the same as in Figure 6. Most of the emission 
measure comes from radii not much larger than Ri^. Therefore, if i?in is much smaller than the 
companion star, most of the emission measure, and all of the emission measure at lower values 
of ^ is occulted during eclipse. The apparent emission measure in eclipse is about one-tenth that 
outside of eclipse here. 

Fig. 8. — The best fit to the ASCA spectra with the stellar wind model. The crosses represent the 
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data. The dashed and dotted lines represent the two absorbed power-law continuum components, 
the Fe Ka line, and the recombination component and the solid line represents the total model 
spectrum. 

Fig. 9. — The eclipse phase model spectrum. The solid line indicates the total flux and the broken 
lines indicate the fluxes of the absorbed power-law continua. Among the features not labeled are 
several Fe 3— >2 lines in the 1-2 keV range. Between Si XIV Ly/3 at 2.375 keV and S XVI Lya at 
2.621 keV is the S XV 2— >1 complex and higher order Lyman lines of Si XIV. 

Fig. 10. — x^ vs. inner wind radius for the disk wind model. The horizontal lines indicate 
confidence intervals. The vertical line represent the neutron star Roche lobe radius, the distance 
from the neutron star to the stellar surface. 

Fig. 11. — Our best fit stellar wind model line fluxes (o) and Ebisawa et al. (1996) measured line 
fluxes (x). Our model line fluxes are always within three sigma of the fit values. All fluxes have 
been corrected for interstellar absorption (A'^hi)- 

Fig. 12. — A comparison of the DEM distributions for Cen X-3 and Vela X-1. 
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